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Competitve hybridization, at the surface and in the bulk, lowers the sensitivity of DNA chips. 
Competitive surface hybridization occurs when different targets can hybridize with the same probe. 
Competitive bulk hybridization takes place when the targets can hybridize with free complementary 
chains in the solution. The effects of competitive hybridization on the thermodynamically attainable 
performance of DNA chips are quantified in terms of the hybridization isotherms of the spots. 
These relate the equilibrium degree of the hybridization to the bulk composition. The hybridization 
isotherm emerges as a Langmuir isotherm modified for electrostatic interactions within the probe 
layer. The sensitivity of the assay in equilibrium is directly related to the slope of the isotherm. 
A simpler description is possible in terms of C50S specifying the bulk composition corresponding to 
50% hybridization at the surface. The effects of competitive hybridization are important for the 
quantitative analysis of DNA chip results especially when used to study point mutations. 

I. INTRODUCTION 

DNA microarrays allow to interrogate the base sequence of DNA or RNA chains. They can be used to detect 
pathogens, identify genetic defects, monitor gene expression etc. (Marshal and Hodgson, 1998; Graves, 1999; Niemeyer 
and Blohm, 1999; Southern et al., 1999; Wang, 2000; Pirrung, 2002). In spite of the intense activity in this field, 
theoretical aspects of the function of DNA microarrays received relatively little attention. Early theoretical work 
focused on the dynamics of hybridization at the surface (Chan et al., 1995; Livshits and Mirzabekov, 1996). Recently, 
theoretical investigations considered the equilibrium hybridization isotherms of DNA chips (Vainrub and Pettitt, 2002; 
Vainrub and Pettitt, 2003) and polyelectrolyte aspects of the systems (Crozier and Stevens, 2003). In the following 
we present a theoretical analysis of the effect of competition between different possible hybridization reactions on 
the sensitivity and specificity of DNA chips. The discussion utilizes hybridization isotherms relating the equilibrium 
fraction of hybridized chains at the surface, x, to the composition of the bulk. The effects are revealed by comparison 
of the hybridization isotherms for competition-free situations with those obtained when competitive hybridization 
is significant. They are quantified in terms of various C50S specifying the bulk composition corresponding to 50% 
hybridization at the surface. A key ingredient of our discussion is the derivation of the competition-free isotherm as 
a Langmuir adsorption isotherm modified to allow for electrostatic interactions. Our model is related to an earlier 
model proposed by Vainrub and Pettitt (VP) in that both assume uniform smearing of the electrical charge of the 
probe layer. 

The elementary units of DNA microarrays are "spots" containing numerous single stranded DNA (ssDNA) chains, 
of identical sequence, terminally anchored to the support surface. The spots are placed in a checkered pattern so 
that each sequence is allocated a unique site. These chains, or probes, preferentially hybridize with free ssDNA 
chains having a complementary sequence. The microarray is immersed in a solution containing labeled ssDNA chains 
whose sequence is not known and are commonly referred to as "targets". The presence of specific sequences is 
signalled by hybridization on the corresponding spot as monitored by correlating the strength of the label signal 
with the position of the spot (Graves, 1999). Recently, label free detection methods, involving optical and mass 
sensitive techniques, attract growing attention (Niemeyer and Blohm, 1999). These allow to monitor the kinetics 
of hybridization. However, such methods measure the total hybridization of a particular probe irrespective of the 
identity of the partner. In marked contrast, selective labeling of a particular sequence monitors only the hybridization 
of this target and does not report on the hybridization of other moieties. 

The unitization of DNA chips as analytical method involves immersing the device in a solution containing a mixture 
of DNA chains of different sequences and concentrations. Under such conditions, it is necessary to allow for the role 
of competitive hybridization. It is useful to distinguish between two types of competitive hybridization. Compet- 
itive surface hybridization occurs when a number of different targets can hybridize with the same probe. Thus, a 
site occupied by certain probes will preferentially hybridize DNA targets with a perfectly matched complementary 
sequence. However, it will also hybridize a certain fraction of mismatched sequences. As we shall discuss, this fraction 
depends on the binding constants as well as the concentrations of the moieties involved. Competitive hybridization 
at the surface clearly lowers both the sensitivity and the specificity of the assay. When the surface competition is 
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FIG. 1: A schematic representation of the competition free case where the probes, p, can hybridize with a single target species, 
t. 
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FIG. 2: In the competitive surface hybridization case the probes, p, can hybridize with a perfectly matched target species, t, 
as well as with a mismatched target, m. 



significant, labeled and unlabeled detection may yield different results. No difference is expected when all targets 
are labeled, as is the case when PCR amplification is used. On the other hand, when selective labeling of specific 
targets is possible, the two techniques measure different quantities corresponding to different isotherms. Competitive 
bulk hybridization reduces the concentration of non-hybridized targets that arc available for hybridization with the 
probe. This takes place when the solution contains complementary sequences that can hybridize with the target in the 
solution. Such sequences may occur either in the same chain, leading to hairpin formation, or in different sequences 
leading to interchain hybridization. Competitive bulk hybridization diminishes thus the sensitivity of DNA chips. Its 
importance varies, again, with the binding constants and the concentrations. The issues discussed above assume their 
clearest form when DNA chips are used to identify single nucleotide polymorphism or point mutations (Lopez-Crapez 
ct ah, 2001). In these situations, the DNA chip is exposed to a mixture of targets differing from each other only in 
the identity of one particular base. The fraction of the different forms is then deduced from the relative intensity of 
the signals of the four spots corresponding to the four possible sequences. 

In practice, the DNA chips are immersed in the target solution for a relatively short time. As a result, the attainment 
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FIG. 3: Competitive bulk hybridization when the probes, p, can hybridize with a single perfectly matched target species, t but 
t can also hybridize in the bulk with a complementary chain, c. c can not hybridize with p. 



of equilibrium is not guaranteed and rates of the different hybridization reactions play an important role. Yet, full 
analysis of the reaction kinetics requires knowledge of the equilibrium state. An understanding of the equilibrium state 
is also necessary in order to identify the relative importance of kinetic and thermodynamic controls of the performance 
of the DNA microarrays. Finally, emerging evidence (Bhanot et al., 2003) suggests that the performance of DNA 
chips, as measured by the number of "false positives" , is best at the thermodynamic equilibrium. With this in mind, 
we investigate the equilibrium hybridization isotherms for three idealized but experimentally attainable situations. 
These situations involve a DNA array immersed in solutions of different composition: (i) A solution containing one 
species of sigle stranded target (Figure 1). (ii) A solution containing two different targets that do not hybridize in the 
bulk but are both capable of hybridizing with the same probe (Figure 2) and (iii) A solution containing two different 
chains, a target and a complementary chain capable of hybridizing with it in the bulk but incapable of hybridizing 
with the probe (Figure 3). In all cases, we consider the case of probes and targets of equal length i.e., the number 
of bases, N, in the chains are identical. For brevity our discussion focuses on systems where the hybridization at the 
surface has a negligible effect on the concentration of targets in the bulk. This case corresponds to small spots or to 
elevated target concentration. 

The first two sections summarize the necessary background information for the subsequent discussion. Thus, section 
II recalls the definitions of sensitivity and other measures of the performance of analytical assays. The relationship 
between sensitivity and the equilibrium hybridization isotherm is also discussed. The relevant structural characteristics 
of DNA chips and important length scales in the problem are summarized in section III. The next section, IV, is 
devoted to the derivation of the competition-free hybridization isotherm as a Langmuir isotherm modified to allow 
for electrostatic interactions. Initially we obtain the hybridization isotherm for an arbitrary electrostatic free energy 
density of the probe layer, 7^. We then consider the hybridization isotherms for particular functional forms of -f e i 
assuming a laterally uniform smearing of the electric charge. We mostly focus on the "diffuse layer" model where the 
charge is uniformly smeared within the probe layer thus allowing for its thickness. It is important to note that some of 
our results are actually independent of the model specifying -f e i . We conclude this section with a discussion of relevant 
experimental results and a comparison between our approach and the VP model. In the remaining sections we pursue 
two complementary goals: The modifications of the hybridization isotherms to allow for competitive hybridization and 
the resulting effects on the sensitivity and specificity of the assay. Three situations are considered. The competition 
free case, when the probes are exposed to a single target, is discussed in section V. This yields upper bounds for 
the sensitivity and the specificity. Competitive surface hybridization is analyzed in section VI and competitive bulk 
hybridization is considered in section VII. The detailed derivation of j e i within the diffuse layer model is described in 
Appendix A. The hybridization isotherm for low salt solutions is discussed in Appendix B. 
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II. ON SENSITIVITY AND THE HYBRIDIZATION ISOTHERM 

As we shall see, the equilibrium hybridization isotherms naturally suggest characterization of the sensitivity of the 
assay in terms of appropriate C50S. This characteristic is closely related to the common definitions of the sensitivity of 
analytical techniques. It is thus useful to first summarize these definitions and their relationship to the hybridization 
isotherms. Different definitions of sensitivity are available (Pardue, 1997; Ekins and Edwards, 1997; references cited 
in Pardue, 1997). The IUPAC definition identifies the sensitivity, S e , with the slope of the calibration curve. The 
calibration curve describes the measured response, R, to a target concentration, ct, R(ct) and 

S e = dR/dct. (1) 

The quantitative resolution of the assay, Act, is then specified by 

Ac t = e r (c t )/S e (c t ) (2) 

where e r is the measurement error as given by its standard deviation. The detection limit, the lowest detectable c*, 
is determined by Ac*(ct = 0) since when the concentration ct is lower than Act(ct = 0) the error is larger than the 
signal. The IFCC convention identifies the sensitivity with the detection limit. 

Our goal is to relate the sensitivity of DNA chips to their hybridization isotherms. With this in mind, it is convenient 
to adopt the IUPAC definition. This choice is motivated by the following observations: (i) the calibration curve in 
equilibrium is closely related to the hybridization isotherm; (ii) The measurement error depends on the measurement 
technique and on instrumental characteristics. In distinction to R(c t ),e r is not related to the calibration curve and 
(iii) <!5 e as given by Eq. H plays a role in the determination of both the qualitative resolution and the detection limit. 

In the following we will assume that R(ct) is proportional to the equilibrium hybridization fraction at the surface, x 
i.e., R{ct) = nx + const where k is a constant. This assumption is justified when the following conditions are fulfilled: 
(i) Non-specific adsorption is negligible and R is due only to hybridization at the surface; (ii) The duration of the 
experiment is sufficiently long to allow the hybridization to reach equilibrium and (iii) the measured signal depends 
linearly on the amount of oligonucleotides at the surface. It is useful to note the following points concerning the 
attainability of these conditions. First, surface treatments repressing non-specific adsorption are available for certain 
substrates (Steel et al., 2000 and references cited in Steel et al., 2000). Second, the attainment of stationary state 
for the hybridization may require long periods of up to 14 hours (Peterson et al., 2001; Peterson et al., 2002; Bhanot 
et al., 2003). Furthermore, the degree of hybridization may depend on the thermal history (heating of the substrate 
or the solution). In this context it is important to stress that, by definition, a state of thermodynamic equilibrium 
is both stationary in time and independent of the path i.e., preparation method. Finally, the linear range varies 
with the measurement technique. For example, when using fluorescent labels the linear regime occurs at low enough 
concentration when self-quenching is negligible (Lakowicz, 1999). 

III. RELEVANT MOLECULAR DIMENSIONS AND LENGTH SCALES 

Two groups of length scales play an important role in our subsequent discussion. One group describes the structural 
features of the probe layer. The second characterizes the electrostatic interactions and their screening. Expression of 
the free energies in terms of these length scales allows for a compact formulation and the identification of the relevant 
dimensionless variables. 

The structural features of the layer are determined mostly by the dimensions of the hybridized and unhybridized 
probes as well as the grafting density (Graves, 1999; Southern et al., 1999; Pirrung, 2002). The number of monomers, 
nucleotides, per probe, N , varies over a wide range. Values of 10 < N < 30 are common but much higher values, 
of N ~ 1000 are attainable. In the following we will consider systems comprising of probes and targets of equal 
size in the range 10 < N < 30. Double stranded DNA (dsDNA) is a semiflexible chain with a persistence length 
ss 10 3 ^4 (Cantor and Shimmell, 1980). Thus, in our N range double stranded oligonucleotides may be viewed as 
rigid rods with the radius of a dsDNA, r — 9.5A and a projected length per monomer along the axis of 26 = 3.4A 
The corresponding parameters for ssDNA are not yet established. Stacking interactions between the hydrophobic 
bases tend to produce a stiff "single stranded helix" (Cantor and Shimmell, 1980; Bloomfield et al., 2000; Korolev 
et al., 1998 and refrences cited in Korolev et al., 1998). Since these interactions are non-cooperative this tendency is 
especially marked in short ssDNA considered by us. Theoretical studies of the melting behavior of free DNA in the 
bulk suggest that ssDNA can be modeled as a rigid rod with projected length per monomer of a « 3AA and a radius 
of r ss « 7 A (Frank-Kamenetskii et al., 1987; Korolev et al., 1998). With this in mind we will approximate the length 
of single stranded chains, Na as identical to that of the double stranded ones, N2b, denoting both by L. For N = 30 
we thus have L 100 A 
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The probes are chemically grafted to the surface via a short spacer chain. The attainable values of the area per 
probe, £, vary with the support surface (Graves, 1999; Southern et al., 1999; Pirrung, 2002). Typical values of E on 
glass surfaces are of order of 10 4 A 2 corresponding to a distance A « 100A between grafting sites. = Significantly higher 
grafting densities of ssDNA are possible on polypropylene supports where E values of E « 40 A 2 , corresponding to 
A « 7 A, were reported. In this last case it is necessary to deplete the surface in order to allow full hybridization to 
take place. The mode of grafting can influence the orientation of the probe. Their orientation can also be affected by 
adsorption to the surface (Levicky et al., 1998). Thus ssDNA grafted onto untreated gold form a compact layer due 
to adsorption. The layer swells and extends into the solution following treatment with mercaptohexanol (Levicky et 
al., 1998). This treatment is also important for elimination of non specific adsorption of the targets. Our discussion 
assumes flexible junctions that enable free rotation and a non-adsorbing surface. Under this conditions, the average 
thickness of the probe layer, H, varies between H w L/2 at low grafting densities and H w L when E<L 2 . 

Three electrostatic length scales are of importance to our discussion. One is the Bjerrum length Is = e 2 /ekT where 
e is the dielectric constant, k is the Boltzman constant and T is the temperature. In water, with e w 80, at room 
temperature, Ib ~ 7 A. Note that the variation of e with T contributes to the T dependence of Is- The second is the 
Gouy-Chapman length A = 1/2-kIb<j. Here a is the number of charges per unit area on a uniformly charged surface. 
A characterizes the spatial distribution of the counterions in the vicinity of a uniformly charged planar surface in a 
salt free solvent. In this situation the majority of counterions are localized within a distance A from the surface. In 
the following the charge of the probes, hybridized or not, is assumed to be uniformly smeared. As a result, a varies 
between N/Ti to 2AT/E depending on x, the degree of hybridization. For an unhybridized layer A is in the range of 
10 — 10 2 A A third scale is the Debeye length, ro, characterizing the screening range of electrostatic interactions in 
a salt solution. For a 1 : 1 salt with number concentration of ions </> s it is m = (87t/b0 s ) -1 / 2 thus, in a IM solution 
r D = 3i. 

The range of DNA concentrations encountered in experiments varies in the range between 10 M to 10 _12 M. The 
solution usually contains also IM of 1 : 1 salt. Under these conditions the electrostatic interactions between the free 
targets are essentially fully screened. 

IV. THE COMPETITION-FREE HYBRIDIZATION ISOTHERM 

The dependence of the hybridization degree, x, on the concentration of the target, c t , is described by the hybridiza- 
tion isotherm. It is helpful to consider first an array of DNA probes of a single sequence, p, in contact with a solution 
containing a single species of ssDNA target, t. The hybridization of p and t creates a double stranded oligonucleotide, 
pt, at the surface. For this choice of system the only reaction is p + 1 ^ pt and no competitive hybridization reactions 
occur (Figure 1). The factors determining the hybridization isotherm fall into two groups. One consists of the factors 
giving rise to the Langmuir isotherm (Evans and Wcnncrstrom, 1994), describing the adsorption of neutral adsorbates 
at a surface comprising a finite number of sites, each capable of accommodating a single adsorbate. These include: 
(i) the entropy of the free targets in solution, (ii) the mixing entropy of the hybridized and unhybridized probes and 
(iii) the non-electrostatic component of the hybridization free energy. The hybridization at the surface of a DNA 
chip differs from the Langmuir scenario in that both the adsorbates (the targets) and the surface (the probe layer) 
are charged. As a result the free energies of the targets and the probe layer incorporate electrostatic terms. These 
allow for the electrostatic interaction energy between the charges and for the entropic effects associated with the 
polarization of the ionic clouds surrounding the macroions. In the following we will obtain a specific form for the 
electrostatic free energy of the probe layer by modeling it as a planar layer with a laterally uniform charge density. 
However, some of our conclusions are actually independent of the functional form of this term. With this in mind we 
introduce at this point an arbitrary electrostatic free energy per unit area 7^. The electric charge localized at the 
surface increases with the fraction of hybridized probes, x. Consequently, j e i = "f e i(x) increases with x, reflecting the 
growth of the electrostatic penalty with the hybridization degree. Initially we will obtain the hybridization isotherm 
in terms of this unspecified 7 e ;(x). We will then consider the hybridization isotherms as obtained for two models for 
the charge distribution within the probe layer and the resulting explicit functional forms of -f e i (x) . 

The equilibrium state of the hybridization reaction, p + t ^ pt, is determined by the condition \i pt = [i v + \i t where 
/ii is the chemical potential of species i. Our discussion focuses on the case where the number concentration of the 
targets is only weakly diminished by this reaction and is well approximated by the initial concentration c t . Since the 
target solution is dilute and the ionic strength of the solution is high, electrostatic interactions between the targets 
are screened. Consequently \i t assumes the weak solution form 

[it = [i° t +kT In c t (3) 

where /j,® is the chemical potential of the reference state. Strictly speaking, \i t = +kT In at where at is the activity 
(Moore, 1972). The dimensionless a t is related to the concentration of t chain ct via at = jet where 7 is the activity 
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coefficient. Since 7 — > 1 as ct — ► we will, for simplicity express [it by Eq. |3| noting that c t in this expression is 
dimensionless. When the concentration of targets is significantly modified by the hybridization with the probes, c t 
should be replaced by c' t = c t — xNt/V where V is the volume of the solution and Nt the total number of probes. 
Such modification is necessary when ct is very low or when the spots are large. 

In order to obtain [i pt we need to first specify the free energy of the probe layer as a function of x. The Nt probes 
are immobilized at the surface thus forming a two dimensional grid of hybridization sites. At equilibrium N p t — xNt 
of the probes are hybridized while N p — (1 — x)Nt remain unhybridizcd. The pt and p chains form thus a two 
dimensional solution associated with a mixing entropy of —kNT[xhix + (1 — x) ln(l — a;)]. This two dimensional 
solution is however non-ideal because of the electrostatic interactions between the chains. Altogether, the free energy 
per probe site is 

Jsite = 70 + x^pt + (1 - x)t4 + ^lei + kT[x\nx + (1 - x) ln(l - x)] (4) 

where E is the area per probe and 70 is the free energy density of the bare surface. n® t and [i® are the chemical 
potentials of the p and pt states in a reference state to be discussed later. For simplicity we now limit the discussion 
to probes and targets with identical number of bases, N. Since each chain carries a charge of —Ne, the number charge 
density on a surface of total area A is a = N(N p + 2N pt )/A — er (l + x) where a = NNt/A is the number charge 
density on the unhybridized surface and E = A/Nt- 

It is convenient to reformulate the equilibrium condition \i pt = fi p -\- (M in terms of the exchange chemical potential 
of the hybridized probe /ipf = fi p t — fi p . The exchange chemical potential of the hybridized probe is \i e p % = dj S it e /dx 
or 

^ =^ t _ fi o + N d^ + kTln _^_ (5) 

where E^r- = E^r^|- = N^^- since da/dx — <tq and Eoo = N. N^^- is thus the electrostatic free energy penalty 
incurred upon hybridization for a given x. The equilibrium condition /xpf = [i t then leads to the adsorption isotherm 
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where Kt = exp ^— A ^ T \ is the equilibrium constant for the hybridization reaction at the surface and AG = 

,,0 _ ,.0 _ 
Ppi rp Pt • 

Our discussion up to this point did not involve a particular model for the charge distribution or a specific functional 
form of 7 e i . In the remainder of this section we will consider the hybridization isotherm for particular forms of "f e i 
as obtained by assuming that the charges of the p and pt chains are uniformly smeared laterally. We will consider 
two models of this type. In the first the charges are distributed in an infinitely thin layer at the solid-liquid interface. 
This model ignores the structure of the probe layer and overestimates 7 e z . It is however of interest as a simple model 
that captures the essential physics. The exact form of 7 e z corresponding to this scenario, for the high salt regime 
encountered experimentally, is specified by the Poisson-Boltzman (PB) equation for rn >C A (Evans and Wennerstrom, 
1994). This 7 e z is identical to the one obtained by the use of the capacitor approximation. In this approximation 7 e ; 
is identified with the electrostatic energy of a planar capacitor, 27r(cre) 2 d/e, with a charge density a = oo(l + x) and 
a width d = rp thus leading to 

^=2^H B r D . (7) 

For this choice of 7 e z the hybridization isotherm Eq. assumes the form 

X =g t exp[-r c (l + a;)] (8) 
ct(l — x) 

where ^7 = T c (l+x) and r c = AirNaolBTD is the electrostatic free energy of a hybridized target in an unhybridized 
layer with a charge density ctq. 

The capacitor model accounts for the essential physics in a simple and transparent way. However this model tends 
to overestimate the electrostatic free energy because all the charges of the DNA chains are placed on a surface. To 
avoid this problem we now assume instead that the charges are uniformly smeared within a layer of thickness H 
giving rise to a number charge density of p = a/H. The analysis of this "diffuse layer" model differs from that of the 
capacitor model only in the form of the electrostatic free energy density 7 e ;. To obtain 7 e ; we utilize a two phase or 
a "box" approximation for the solution of the PB equation (Pincus, 1991; Wittmer and Joanny, 1993; Borisov et al., 
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1994). Within it, we distinguish between two regions: (i) a proximal region, adjacent to the charged surface, where 
the concentrations of ions deviates from the bulk values. The concentrations of each of the ionic species are constant 
and obey the Donnan equilibrium, (ii) A distal neutral region where the effect of the charged surface is screened out 
and the concentrations of the ions are determined by the concentration of the salt. The ionic concentrations and the 
equilibrium electrostatic free energy are determined by minimization of the free energy with respect to the height of 
the proximal region. This approximation involves the simplest form of discretization of the PB equation. The details 
of the analysis are presented in the Appendix A. In the following we focus on the experimentally relevant case of high 
salt such that td <C H and <C (HA) 1 / 2 . The low salt regime is described in Appendix B. In the high salt regime 
the screening of the charged layer is dominated by the contribution of the salt and 

l± = ^l B r \. (9) 
kT H y ' 

The hybridization isotherm in this "salt screening" (ss) regime is 

Ct{1 X _ x) =^cxp[-r(i + x)] (io) 

where -^p^gf- ~ SnNals-ff- — T(l + x) and V = SttNo-oIb^j- is the electrostatic penalty incurred by a pt chain 
in an unhybridized layer with a = <tq . Note that the functional form of Eq. ^| is identical to that of Eq. [H] but 

r = 2T c r D /H < r c . 

As a reference state it is convenient to choose the state of a chain (ssDNA or dsDNA) anchored to a surface at a 
low grafting density such that the in-plane electrostatic interaction are negligible. When the "lateral" interactions are 
negligible, one may roughly approximate jJl t (/4J) by the fjP of the corresponding free chain in the solution. This choice 
is useful in that it enables us to estimate the various hybridization constants using the nearest neighbor parameter sets 
available in the literature (Bloomfield et al., 2000). It is however important to keep in mind the problems introduced 
by this choice of reference state and the approximation of /j,® t (/i°). One difficulty involves the electrostatic free energy. 
7ei is obtained by charging of a hypothetical non charged layer. As a result, the electrostatic contribution to fi® t 
leads to a small overestimate of the electrostatic free energy. Note that for high a or small A fluctuation effects become 
important (Lau et al., 2002). These are not included in our analysis. In addition, caution is required in identifying 
the boundaries of the regime of negligible lateral interactions. This is because the decay of electrostatic interactions 
at an impenetrable surface is slower than in the bulk. Thus, point charges embedded at an impenetrable surface 
polarize an hemisphere of the ionic solution thus giving rise to a dipole and the lateral interactions decay as 1/r 3 
(Jancovici, 1982). Another problem concerns the rotational free energy of the chains. The rotational freedom of the 
terminally anchored chains is restricted by the impenetrable grafting surface. Further restrictions may be imposed by 
the grafting functionality. The diminished rotational freedom reduces the rotational term in the free energy per chain. 
This effect is however neglected when \j^ t (/z°) are approximated by /jP of the corresponding free chains. When both 
the target and probe are self complementary it is necessary to allow for the change of symmetry due to the grafting. 
In turn, this requires an appropriate modification of /i° t (/x°) with respect to their bulk counterparts. Finally, note 
that in the low grafting density regime, as discussed above, the hybridization isotherm is expected to assume the 
Langmuir form 

K t . (11) 



(1 - x)ct 



In this regime the electrostatic aspect of the problem is evident only in the dependence of the /j,°s and thus K t , on 
the concentration of salt. 

The number of hybridization isotherms of DNA chips reported in the literature is rather small (Nelson et al., 
2001; Peterson et al., 2001; Peterson et al., 2002). The situation is further complicated because of paucity of data 
concerning Nt, the number of probes available to hybridization, and the related problem of ascertaining the attainment 
of thermodynamic equilibrium. The "uniform smearing" models for the hybridization isotherms are supported by two 
experimental studies carried out by the group of Georgiadis. In one experiment the grafting density was varied in 
the range of 2 • 10 12 — 12 • 10 12 probes/cm 2 while Ct was kept constant at 1/iM (Peterson et al., 2001). A plot of 
ln[(l — x)/x] vs. (1 + x)/Tj can be fitted with a straight line with a slope smaller than the one predicted by the 
theory (Figure 4). This is however encouraging since the data was acquired in l/2h and is thus unlikely to reflect 
complete equilibrium. In the second group of experiments, the hybridization was studied for a lower grafting density 
of 1.5 • 10 12 probes/cm 2 while Ct was varied over the range of 500mM to 5/xM (Peterson et al., 2002). In this study 
the hybridization isotherm of the perfectly matched targets was well fitted by the Langmuir form. Importantly, this 
study established that the system failed to reach equilibrium without heating treatment for mismatched targets. 
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FIG. 4: Aplot ofln(l-as)/a: vs. (l+z)/£ using the data of Peterson et al. (2002). Ea.Hulvields InCl-xl/x = const' +B(l+x) IT, 
with B = 8Trl B N 2 r 2 D /H. For the experiment cited l B = 7 A, r D = 3A, N = 25 and H = 85i leading to B ~ 1.16 x 10 4 i 2 as 
compared to the observed B~3x 10 3 ^4 2 . 



A hybridization isotherm of identical form to Eq. |H1 and to Eq. 1101 was announced earlier by Vainrub and Pettitt 
(Vanirub and Pettitt, 2002; Vanirub and Pettitt, 2003). VP also pointed out that some of the results of the Georgiadis 
group are consistent with this form. The VP approach is designed to permit the utilization of exact results on the 
interaction free energy between a penetrable charged sphere and an impenetrable charged surface in the strong 
screening regime when the Debye-Hiickel approximation is applicable (Ohshima and Kondo; 1993). Within it, one 
calculates the excess free energy of a probe layer with xNt hybridized probes, F e i(x), with respect to the unhybridized 
layer. In effect, F e i(x) is the sum of the contributions of xNt hybridization events, F e i — -^iO 7 *)- Each step 

contributes Fi(ai) = F pt (<Ji) — F p (<Ti) where F pt (F p ) is the electrostatic free energy of a pt (p) sphere in contact 
with a planar layer with charge density o~i = o"o + iN/A. Thus, at each step the probe layer is modeled as a planar 
charged surface interacting with a single charged sphere. The steps differ in the charge density of the surface. The 
main difference between the VP approach and ours is in the handeling of the charges. In the VP scheme some of the 
charges appear as charged spheres while others appear as a charged surface. Within our model there is no duality and 
all charges are described in the same fashion. In practical terms, the VP approach can not allow for the thickness of 
the probe layer nor can it be extended to describe hybridization at lower ionic strength. 



V. SENSITIVITY, SELECTIVITY AND cso FOR COMPETITION FREE SYSTEMS 

The hybridization isotherms discussed in the two preceding sections describe DNA arrays in the absence of compet- 
itive hybridization in the bulk or at the surface. This situation is realized when an array comprising of single type of 
probes is exposed to a solution of a single target. The concentration of target leading to 50% equilibrium hybridization 
in such systems, 'cgp = K^ 1 exp(^^ i | :c=1 /2) is a useful characteristic of the system. Within the diffuse layer model 
in the salt screening (ss) regime t c^ >Q is 

t C ° = -^exp(|r). (12) 



4 is closely related to the sensitivity of the array 



s *> " i + < I("i-l ) r K " B " > '- r " + I " " i + ",("i-'.)r 4" p '" r(x " 1/a)1 - <13) 

The sensitivity of the array, as defined by S e {x) varies with x and thus with ct- It is maximal at x — when 

S e (0) = fftexp(-r) = r i r exp(-^) (14) 

c 50 z 
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while at x = 1/2 it is 56(1/2) = 1/(4 + r)*c§ . As we shall see, 5 e (0) is not affected by competitive hybridization. 
On the other hand, S e (x) and C50 are modified significantly by these processes. 

Since S e (x) ~ 1/*C5 , clearly a lower t c^ is desirable and l/'cgg is a useful measure of the sensitivity of the array. 
Both l/*Cg and S e (x) decrease as T and the electrostatic penalty incurred by the hybridization increase. In the salt 
screening regime, where most experiments are carried out, T increases with the grafting density as T ~ <jq. While 
higher sensitivity is expected at lower grafting densities, this does not ensure a lower detection limit or a better 
quantitative resolution. These last two parameters depend also on the measurement error e r . In turn, e r typically 
decreases as the grafting density, and the signal, increase. Thus, l/*Cg and S e (x) only provide partial guidance for the 
design of DNA arrays. Nevertheless, these two parameters do provide useful information regarding the performance 
of a DNA chip of a given design (that is, grafting density, grafting functionality, spot size and detection method). 
Thus, the relative sensitivity of two different probe target pairs, p\ t\ and P2 t2, all other factors being equal, is 



L o t 2 „o TS 



*2 Q H „U js 
D e C 50 J\t 2 



(15) 



The specificity of a given probe, p, can be quantified by the relative sensitivity when a p spot is exposed to a perfectly 
matched target, t, or to a mismatch, m 



mj) 



(16) 



m Q t„0 IS " 

These two ratios also specify the corresponding ratios of the qualitative resolution and the detection limit. Importantly, 
Eq. E| and Eq. ^| are independent of the electrostatic penalty irrespective of the form of -f e i . 

VI. THE EFFECT OF COMPETITIVE SURFACE HYBRIDIZATION 

The hybridization isotherm requires modification when the bulk solution contains more than one ssDNA species 
capable of hybridization at the surface. In this situation the different species compete for hybridization with the 
probes. For simplicity we consider the case of a binary solution comprising a target (t) and a mismatched ssDNA (m) 
with a concentration c m and a standard chemical potential in the bulk solution /x m . It is placed in contact with a 
single component probe layer such that the p chains are perfect matches to the targets (Figure 2). We further assume 
that the m and t chains are of the same length. The number of probes that hybridized with m is N m — yNr- In this 
case a = N(N P + 2N pt + 2N pm )/A = ct (1 + x + y) and 

Isite = 7o + x^ip t + y^p. m + (l - x - y)^,p + £7^ + kT[x\nx + y lny + (1 - x - y) ln(l -x-y)} (17) 

where yP pm is the standard chemical potential of a hybridized pm at the surface. The hybridization isotherm in this 
situation is determined by two equilibrium conditions = fit, as before, and \i e v x m = fx m . In obtaining the explicit 

form of these conditions note that d ^p x el — d ^Qy l = N^^- because |~ = ^ = <7q. The exchange chemical potentials 
of the hybridized m and t are thus given by 

& - & ~ M° + + fcTln T -^— (18) 

^ = $m -nl + N^ + kT\n l _ V x _ y (19) 
and the chemical potential of the free m is 

Mm = Mm + kThiCm. (20) 

As before, we focus on the "small spot" limit where the bulk concentrations of m and t are not affected by the 
hybridization at the surface. The hybridization behavior of this system is described by three isotherms specifying the 
hybridization degrees of t and m individually as well as the total hybridization: 

X = Ktexp " N d lel / 21 s 

ct(l — x — y) ' kT da J 
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FIG. 5: Plots of 'csoAso vs. c m / 



TO o 
c 50i 



as given by Eq. 1251 for the case of competitive surface hybridization involving 



a probe, p, of the sequence CAACTTGATATTAATA, a target, t, GTTGAACTATAATTAT and a mismatched target, m, 
GTTGAGCTATAATTAT (TG mismatch). In the three cases depicted T = 300°^, N = 16, H = MA, l B = 7A and r D = 3A. 
The continuous line corresponds to the low grafting density regime where T = 0. The two other are E = H 2 = 2916i 2 leading 
to r = 2.57 (dashes) and E = lO 3 ^ 2 leading to F = 7.5 (dots). The standard Gibbs free energies per mole at 37° G are 
AG? = 12Akcal/mole and AG™ = 10. lkcal/mole (Tibanyenda et al., 1984). Since the AG are per mole rather than per 
molecule, the equilibrium constants at T — 300° K, neglecting the T dependence of the AG , are Kt = exp(— AGt/RT) ~ 10 9 
and K m = exp(-AG„/i?T) ~ 10 7 ' 4 where R is the gas constant. The corresponding t c% are 10" 9 M, 10~ 7 ' 4 M and 10" 4 A M 
respectively. The values of m c% are 10" 7 ' 4 M, 10" 5 7 M and 10" 2 ' 5 M 



y 



c m (l -x-y) 



K m exp 



N d le i 
kT da 



(22) 



{I- x-y) 



= (c m K m + ctK t ) exp 



N d le 



kT da 



(23) 



where K t — exp {^—^fr-^j , K m — exp A ^"' j and AG° ra = ^ pm — /i° — $ n . The observed isotherm depends on the 

method used to interrogate the surface. Thus, utilization of selectively tagged t will reveal Eq. use of selectively 
tagged m will show Eg. 1221 while detection methods sensitive to overall hybridization mass, such as surface plasmon 
resonance, will yield Eq. 1231 The explicit form of the hybridization isotherms within the diffuse model in the salt 
screening regime is obtained by substituting ^7^- = r(l + ir + y). Note that K t , K m and T can be determined from 
experiments involving exposure of the DNA chip to single component solutions of t and m chains. 

The specificity of the assay can be quantified by considering the fraction of "incorrectly" hybridized probes, P m . 
Equations 1211 and 1221 yield y = x — and thus 



P„ 



y 



c t K t 



(24) 



Within this definition the specificity strongly depends on c m , or to be precise on 



probes is small, P m < 1, so long as c m <C c t 



K,„ 



Ct 



c m K„, 
c t K t 



The fraction of mismatched 



, half of the hybridized probes are mismatched, 



Pm = 1/2, while for c rn Ct^jP^, P m approaches unity. Equation 1241 is independent of the electrostatic contribution 
irrespective of the form of 7^. It is also useful to consider the ratio of t c\ Q to '050, the bulk concentration of t giving 
rise to 50% pt hybridization in the presence of a mismatch of concentration c m . In contrast to P m , the expression for 



1 Ji 



c 5o/* c 50 does depend on 7^. For the diffuse layer model in the salt screening regime it is given by 
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In the low grafting density regime, when T = 0, Eg. 1251 assumes the form -ffi- = 1+ -^rr- In all cases, *cso = *Cr 

when c m — and 'cso > 'cg for c m > 0. In other words, the sensitivity, as measured by l/'cso, decreases as c m 
increases (Figure 5). 

VII. THE EFFECT OF COMPETITIVE BULK HYBRIDIZATION 

A different type of competition occurs when the targets can hybridize in the bulk as well as at the surface. Such 
competition can arise in three different situations: (i) The solution contains targets as well as complementary strands, 
c. These can be perfectly matched or mismatched sequences. The c chains hybridize with the targets to form free 
double stranded tc DNA chains. Thus, the t + c ^ tc reaction in the bulk competes with the t + p ^ pt reaction at 
the surface (Figure 3). (ii) The targets are self-complementary and thus capable of undergoing a bulk hybridization 
reaction t + t ^ tt in addition to t + p ^ pt where p now denotes the immobilized t probe, (iii) A third possible 
scenario involves formation of hairpins. 

As explained in section III, within our discussion the lengths of the p and pt chains are identical. Accordingly we 
will focus on the first two cases where the length of the chains does not change upon hybridization. Initially, we 
discuss the t + c ^ tc scenario and then comment on the modification required to adapt the analysis to the t + t ^ tt 
case. Again, we focus on the "small spot" limit assuming that the hybridization with the probes has a negligible 
effect on the concentration of the targets. The hybridization isotherm describing this situation, for the two cases of 
interest, is 



(1 



x is 
— = K t exp 



N d lel 
kT da 



(26) 



and -^^g^ = T(l + x) in the ss regime of the diffuse layer model. Importantly, the hybridization isotherm is modified 
in that q, the total concentration of t, is replaced by the equilibrium t concentration, [t]. In turn, [t] is determined 
by the mass action law governing the bulk hybridization reaction. The combination of Eg. 1261 with the appropriate 
mass action law is equivalent to the equilibrium condition specified by fit + fJ-p — Upt and fj, t + fi c — fJ-tc- 

In the t + c f± tc scenario the mass action law is [fc]/[i][c] = K where [i] is the equilibrium concentration of 
species i and K is the equilibrium constant of the bulk hybridization reaction for the temperature and ionic strength 
considered. This is supplemented by the mass conservations relations [t] + [tc] = Ct and [c] + [tc] = c c where Ci denotes 
the total concentration of i. [t] is then specified by 

K[t] 2 +{K(c c -ct) + l}[t]-ct = 0. (27) 

When the hybridization with the probes has a significant effect on the concentration of the targets, [t] + [tc] = c t 
should be replaced by [t] + [ct] + xNt/V — c t . For brevity, we will not consider this case. It is instructive to analyze 
the effect of the competitive bulk hybridization for a number of simple situations. When the equilibrium favors the 
reactants, [t] w c t and the hybridization isotherm retains the competition-free form, Eq. Such is the case in the 
presence of large excess of t, ct 3> c c or when K is sufficiently small i.e., c c Ct or c c rj c* but Kc c ^ 1. Significant 
modification of the hybridization isotherm occur when the bulk hybridization equilibrium favors the products. This 
situation occurs in two simple cases: when Kc c 3> 1 with either c c ^ Ct or c c ~ Cf. We initially discuss briefly the 
first situation when 

M « Yc~ « c * (28) 



leading to 



(1 - x) Kc, 



K t exp 



N d le 



kT da 



(29) 



To obtain an explicit form of the isotherm within the ss regime of the diffuse layer model we substitute T(l -I- x) for 



j^^O^-- However, the effect on 'cso is independent of the model. In comparison to t c^ = K t 1 exp ^-^^^-U^ 



1/2 



*C5o increases to 



{ c 50 = Kc^o » s ^ . (30) 

The sensitivity, as measured by l/*C5o, is thus reduced by a factor of A'cc ;§> 1. When c c ~ ct and A'cc ^> 1 the 
equilibrium condition H27(l yields 

•C t \l/2 
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thus leading to 



\k) Ktexp 



(1-x) 

The corresponding *C5o increases thus to 



n d lel 

kT da 



(32) 



c 50 = KCcl ) 2 (33) 



and the sensitivity, as measured by 1/*C50, is reduced by a factor of -f^'c^Q 3> 1 in comparison to the competition-free 
scenario. The sensitivity S e — dx/dc t does depend on the form of j e i. When Eg. 1331 is applicable S e , as specified by 
the uniform density model at the ss regime, is 

K 2 (1 — x) 3 1 (1 — x) 3 

S e = -^- exp[-2r(l + a;)] M V = ./ rr = ; exp[-r(x - 1)1 .„ y _ > -r. (34) 

AK L v n x[l+Tx(l-x)] *c 50 n x[l + Tx(l-x)} y ' 

However, in the limit of c± — ► the effect of the competitive bulk hybridization is negligible and i5 e (0) is thus given 
by Eq. [21 This is also the case for the c c 3> c% and K c c ^> 1 scenarios considered earlier. 

In the low grafting density regime, when j e i is independent of a, the hybridization isotherm for c c w c t with 
Kc c 3> 1 assumes the form x/(l — x) = K t (ct/K) 1 / 2 . Upon defining if e // = K 2 /K this isotherm can be expressed 
as 

« . (35) 

This form is of interest because it resembles the isotherm obtained from the Sips model (Sips, 1948). The Sips model 
provides a generalization of the Langmuir isotherm in which the single binding energy, utilized in the Langmuir 
version, is replaced by a distribution of binding energies thus leading to an expression of the form 

- {Ke " Ct)a (36) 



1 + (K eff c t Y 

where a is a characteristic of the distribution function. Thus, competitive bulk hybridization can give rise to a "Sips 
isotherm" with a = 1/2 even though the underlying mechanism is completely different. This is of interest because the 
Sips isotherm was recently reported to allow for improved fitting of hybridization data (Peterson et al., 2002). 

When the competitive bulk hybridization involves self-complementary chains, t + t f± tt, the preceding discussion 
requires modification. In this case the mass action law assumes the form [tt]/[i] 2 = K and the corresponding mass 
conservation relation becomes [t] + 2[tt] — c t . [t] is thus determined by 2K[t] 2 + [t] — c t — 0. When Kc t <§C 1 the 
competitive effect is negligible and [t] ~ Ct- In the opposite limit, Kc t 3> 1, the bulk hybridization is important and 
[t] fa {ctjlKfl 2 . The t + t ^ tt scenario thus closely resembles the t + c f± tc case when c t pa c c . Note however 
that care must be taken in estimating K t for the self complementary case. When the sequences of the p and t chains 
are identical, K t differs from the bulk K because the grafting to the surface modifies the symmetry of the chain (in 
addition to the factors discussed in section IV). 



VIII. DISCUSSION 



The hybridization isotherms of DNA chips provide a natural starting point for the analysis of their sensitivity and 
specificity. Clearly, this description is incomplete in that it is limited to equilibrium states while in typical experiments 
equilibrium is not attained. The hybridization isotherms are nevertheless of interest because of the emerging evidence 
that the best performance of DNA chips is obtained in thermodynamic equilibrium (Bhanot et al., 2003). Accordingly, 
the selectivity and specificity obtained from the hybridization isotherms provide upper bounds to the performance of 
these assays. This approach is also of interest because an understanding of the equilibrium state is a prerequisite for 
the full analysis of the kinetics of hybridization. When selectivity is discussed in terms of the slope of the response 
curve, it is necessary to use an explicit form of the hybridization isotherm. We obtained such an explicit expression 
by use of the diffuse layer model. In this model the charges of the pt and p chains are uniformly smeared within the 
probe layer. However, the analysis of the hybridization isotherm also suggests the use of various C50S as measures of 
the specificity and selectivity of DNA chips. This description affords an important advantage in that the effects of 
competitive hybridization can be described in a form that is independent of the model used to specify the electrostatic 
interactions. Thus, the best performance of DNA chips is attained in competition-free situations used to define 'c^q, 
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m Cg etc. One can then analyze the effects of competitive hybridization in terms of the increase in 'cso in comparison to 
*Cg . This analysis also indicates that the knowledge of the competition- free isotherms allows to predict the isotherms 
realized when competitive hybridization occurs. In addition the observed isotherm depends on the measurement 
technique when competitive surface hybridization is important i.e., label free detection differs from the detection of 
selectively labeled targets. 

Much of our discussion concerns the effects of competitive hybridization. In certain applications the effects of 
competitive surface hybridization can be minimized by proper design of the probes (Lockhart ct al., 1996; Li and 
Stormo, 2001; Bhanot et al., 2003). Such is the case, for example, when studying the expression level of genes of 
known sequence. However, this strategy can not be employed when DNA chips are used to identify single nucleotide 
polymorphism or point mutations. Probe design is also of limited value in counteracting the effects of competitive 
bulk hybridization. 

The results we obtained are based on the equilibrium hybridization isotherms. They are formulated in terms of 
the equilibrium fractions x, y etc. of hybridized probes. In confronting these predictions with experimental results 
it is important to note the following two points. First, in order to specify x and y it is necessary to determine first 
the number of probes available to hybridization, Nt- Thus it is not sufficient to ascertain the number of p chains 
immobilized at the surface. It is also necessary to confirm that this corresponds to the number of hybridized probes 
at equilibrium with a large excess of targets. This brings us to the second point concerning the equilibrium state. 
This plays a role both in the determination of Nt, as discussed above, and in the determination of equilibrium 
fractions of hybridized probes. Here we recall again that a stationary state does not necessarily imply equilibrium. 
An equilibrium state should also be independent of the preparation method or sample history. In the context of DNA 
chips it is thus important to verify that the stationary state is not affected by a heating treatment. In every case, 
the equilibration time can be very long with periods of up to 14 hours reported in the literature. It is also useful to 
note that the equilibration time depends on the bulk composition, ct and c m , on the ionic strength and the grafting 
density, S. It also varies with the number of mismatches and their identity. Accordingly, the equilibration time in one 
experimental situation is not necessarily identical to the equilibration time under different conditions. When studying 
simultaneously the hybridization on different spots the equilibration rates for the different spots may well differ. 

It is useful to distinguish between two types of experiments involving DNA chips: experiments designed to elucidate 
the physical chemistry of their function and experiments utilizing DNA chips to analyse biological samples. In the first 
category, the experimental set up allows for selective labeling and for precise control of the composition of the bulk 
solution. It is straightforward to confront our analysis with such "physical chemistry" experiments. The situation 
with respect to analytical applications is more complex. Analytical experiments typically rely on PCR amplification of 
biological samples. As a result, selective labeling is impossible and the composition of the bulk solution is determined 
by the composition of the original sample and the amplification scheme i.e., the choice of primers. Our discussion 
reveals difficulties in the quantitative interpretation of the results of such experiments, especially when used to study 
point mutations. In this last situation, one may quantify errors introduced by the competitve hybridization by use 
of "standard addition" i.e., study a series of solutions obtained from the amplified biological sample by addition 
of different amounts of synthetic, selectively labeled target. The practical importance of these difficulties and the 
methods to overcome them remain to be established. 



IX. APPENDIX A: THE BOX MODEL FOR A DIFFUSE AND FOR A PLANAR LAYER. 

We consider a diffuse layer carrying Q charges distributed uniformly in a region of height H such that the total 
charge is — Qe < 0. The resulting number charge density is p = Q/AH = a/H where a = Q/A is the corresponding 
surface number density of charges and A is total surface area. In the limit of H — this system reduces to the case 
of a charged surface. The analytical solution of the PB equation for this last case is known. Accordingly we will also 
investigate the box model for the H = in order to demonstrate that it recovers the known results up to a numerical 
factor. 

The surface charge affects the distribution of ions within a proximal layer of height A > H , adjacent to the surface. 
Within this layer n± is the total number of univalent positive (negative) ions and <j>± = n±/\A are the corresponding 
number concentrations. The electrical potential in the "box" , \&, determines the deviation of <j>± from the bulk number 
concentration <p s via <j>± — 4> s exp(±e v I' / kT) thus leading to the Donnan equilibrium 




2 



(37) 



8 ' 



The overall electroneutrality of the proximal layer, n + — n_ = Q leads to 



A(j) = <p + - <p- = a/X. 



(38) 



14 




X 



FIG. 6: The concentration profiles of ions within the box model for the diffuse layer. The uniformly smeared charge of the p 
and pt chains is depicted by the shaded step function. It causes the concentration of negative and positive ions, <f>- and <j>+, 
within the proximal layer of thickness A, to deviate from the bulk value 4> s . 



A is the "neutralization length" of the system in that the net charge of a thicker layer is zero and at higher altitude 
vjr = o. Combining Eq. l3"7l and Eq. |2H1 leads to a quadratic equation <j^_ — j(f) + — <j) s = determining <j>±. Upon 
introducing the parameters s = rp/A and x = A/ A we obtain 



±- 



2.s 2 



4s 4 



1/2' 



(39) 



The excess entropy of the ions in the box, with respect to the bulk, is specified by — S/k = n~ ln(^>_/<^ s )+n+ ln(0+ /4> s )- 
Invoking Eq. 1371 and Eq. 1381 leads to —S/k = Aa\n{<j) + / <j) s ) and the excess entropy per unit area is thus 



S_ 

Ak 



a In 



2s* 
x 



4s 4 



1/2' 



(40) 



The charge per unit area that is bound by a surface of height z is ez(— p + A</>) when < z < H and e(— p + zAcf>) 
when H < z < A (Figure 6). Consequently, the electrostatic field, E(z), as determined by the Gauss theorem, is 



E{z) 



E in {z) = 
E out {z) = 22* 



("1 + f 



< z < H 
H < z< A 



(41) 



In the H = case the charge per unit area below z is e(— p + zA(j>) and E(z) = E ou t(z) for < z < A. The associated 
electrostatic energy per unit area, W = ^ J E 2 (z)dz is 
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In the case of H = this reduces to ax/3. Altogether, the electrostatic free energy per unit area is 



7 ei ax f H 

fcr = irl iea + (Tln 



The equilibrium condition d"/ e i/dx = leads to 



1 + 



4s 4 



4\ V2 



We first consider the H — case when 



z z 1 



4s 4 



2s 2 f As" 

1 " 



H 

xA 



4\ 1/2 



= 6s 2 



4\ 1/2 



6s 2 . 



(43) 



(44) 



(45) 



In the high salt limit, when s <C 1 this leads to equilibrium values of x w 6 1 / 2 s and j e i/kT w 2(2/3) 1 / 2 crs « 1.6<ts or 



A«6 1/2 r D 7e; //cT«47r(2/3) 1/2 CT 2 ; 



(46) 



as compared to ^ e i/kT — as obtained from the rigorous solution of the PB equation. In the opposite limit, of s 3> 1 
corresponding to low salt, Eq.QSlleads to x » 3 and j e i/kT « 2er[ln2s + (1 — ln3)/2] or 



A w 3A 7 ei //cT « 2cr ln(47ro-Z B r£)) 



(47) 



while the rigorous solution of the PB equation is jei/kT « 2cr[ln2s — 1]. Thus, the box model for the planar layer 
recovers the rigorous solutions of the PB equation up to numerical corrections. In the low salt regime it yields 
the correct leading term j^/kT ss 2c Ins. However, at high salt the box model overestimates j e i by 60%. This 
performance is indicative of the errors expected from the model for the diffuse layer. 

When H > the equilibrium condition Eq. ^]is applicable. This equation differs from Eq. 03] in two respects: 
(i) a [l — (H/xA) 2 ~\ factor arising from the modification of the charge distribution and the associated electrostatic 
energy and (ii) the problem now contains an additional length scale, H. We expect that A > H and consequently 
the magnitude of 4s 4 /x 2 = 4r|j/A 2 A 2 can be large (small) even when s = t\d/A <C 1 (s 3> 1) provided U < rjj 
(H 3> re,). To allow for this last feature it is convenient to express Eq. ^]in terms of y = X/H instead of x leading 
to 



(y 2 i) 



4s 4 

y 



1/2 



= - > £ 



(48) 



In analyzing the asymptotic solutions of this equation it is useful to compare the neutralization length, A, with 
H. Two principle regimes emerge. When A > > 1), the structure of the diffuse layer is irrelevant and we 

recover the solutions of the PB equation describing a charged planar layer. In this "PB limit" Eq. 021 reduces to 

6s 2 (A) 2 . Here we can again distinguish between two regimes. When s 2 A/yH ^> 1 this leads 



y 



l 



r(£) S 



1/2 



to y m 3A/H > 1 while for s 2 A/yH «lwe obtain y rj 6 1/2 sA/H. Altogether 

A; 



3A 

6 1 / 2 riD 



H < A and r D > A 
H < and r D < A 



(49) 



When A ks A the screening of the electrostatic potential is due to the counterions of the charged layer. The coions, 
originating from the salt, dominate the screening when A ~ rjj. The crossover between the "salt screening" (PBss) 
and "counterions screening" (PBcs) regimes in the PB limit occurs at s 2 A/yH = 1 leading to s — 1 or A = rjj. 

When y > 1 the charge distribution within the diffuse layer plays an important role. In this case it is useful to 
express y as y = 1 + 6 and to solve with respect to 5 < 1. Eq . EEl reduces to 2 8[1 + (2 s 2 A/H) 2 } 1/2 = 6s 2 A 2 /H 2 . 
Consequently we can distinguish between two cases depending on the magnitude of s 2 A/H. When s 2 A/H ^> 1 or 
r 2 D > AH we obtain S ps 3A/2H. In the opposite limit, of s 2 A/H < 1 or r 2 D < AH we obtain 5 w 3 {A/H) 2 s 2 . That 



is 



H + §A H > A and r 2 D > AH 
H + 3^ H^>r D and r 2 D < AH 



(50) 
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FIG. 7: The asymptotic regimes of the diffuse layer within the box model. In the two PB regimes (PBcs and PBss) the 
neutralization length A is large, A 2> H, and the layer behaves as a charged planar surface. In the two remaining regimes, 
A > H and the charge distribution of the layer, p, plays a role. In the cs regions the screening is dominated by the counterions 
while in the ss regions it is due to coions originating from the salt. 



When A w H + 3A/2 the screening is due to the counterions while for A m H + 3rjj/H it is dominated by the coions. 
The crossover between the "salt screening" (ss) and "counterions screening" (cs) regimes is specified by s 2 A/H = 1 
or A = r 2 D /H. Additional crossovers clearly occurs at rrj = H and at A = H (Figure 7). 

To obtain the corresponding asymptotic expressions for 7^ it is convenient to rewrite Eq. 1481 in terms of y as 

1/2" 
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kT 



H (y - 1) 
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yH 



yH 



(51) 



When A/H <C 1 and s 2 A/H 3> 1 (cs regime), y as 1 + 3A/2H, the logarithmic term is dominant and j e i/kT « 
crln(4s 2 A/i7). In the limit of r D /H <C 1 and s 2 A/H <C 1 (ss regime), when y w 1 + Sr^/H 2 , the logarithmic term 
can be expanded in powers of s 2 A/H leading to j e i/kT ~ 2ar 2 D /HA. When A/H 3> 1 and s 2 A/yH ^> 1 (PBcs 
regime), the logarithmic term is dominant and jei/kT « er[l + ln(4s 2 /3)] while for ro/H 3> 1 and s 2 A/yH <C 1 
(PBss regime), the logarithm can be expanded leading to i e i/kT « 2 (2/3) 1 / 2 crro /A. The four scaling regimes are 
summarized in Table I. 



cs 


kT 

aln (al B r 2 D /H) 


H + 3A/2 


A < H and r D > {AH) 1 ' 2 


ss 


<? 2 l B r 2 D /H 


H + 'ir 2 D /H 


r D < H and r D < {AH) 1 / 2 


PBcs 


a\ii{(jl B rD) 


3A 


A > H and r D > A 


PBss 


a 2 l B r D 


6 1/2 r D 


r B > H and rjj < A 



X. APPENDIX B: THE HYBRIDIZATION ISOTHERM AT LOW SALT. 



A novel form of the hybridization isotherm is obtained at low salt, when the screening is dominated by the coun- 
terions of the p and pt chains. This is the case when the concentration of counterions within the probe layer is much 
larger than the concentration of coions contributed by the salt leading to > (AH) 1 / 2 and A ^ H. In this situation 



kT V H 
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The hybridization isotherm in this "counterion screening" (cs) regime is 



x 



= K t exp [-r„-JVIn(l + a;)] 



(53) 



ct(l - x) 




where ~ r cs + ATn(l + x) and T cs = N[\n [87ra lB-ff-J + 1]- The cs regime is of interest in that it provides 

an additional test for the diffuse layer model. 

The authors benefitted from instructive discussions with T. Livache and P. Pincus. EBZ was funded by the CNRS 
and the Universite Joseph Fourier. 
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